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ABSTRACT 


A  new  nonlinear  time-domain  method  for  predicting  large  amplitude 
motions  for  a  ship  advancing  in  a  seaway  has  been  developed.  The  seaway 
is  represented  by  a  second-order  Stokes  wave.  A  consistent  nonlinear  theory 
is  derived  by  assuming  that  the  frequency  of  encounter  is  small.  In  this 
theory,  the  added-mass,  damping,  and  diffraction  terms  are  obtained  by  linear 
theory,  whereas  the  nonlinear  hydrostatic  restoring  and  Froude-Kriloff  forces 
are  predicted  by  integrating  the  pressures  over  the  instantaneous  wetted  sur¬ 
face  of  the  hull.  Computer  codes  have  been  developed  which  predict  the  non¬ 
linear  large-amplitude  motions  in  the  time  domain  for  a  ship  advancing  in 
beam,  head,  and  following  seas.  Results  obtained  by  these  computer  codes  show 
good  agreement  with  linear  theory  for  small -amplitude  waves,  whereas  large 
differences  are  found  between  the  nonlinear  method  and  linear  theory  for 
steeper  waves. 


Section  1 


INTRODUCTION  AND  SUMMARY 


1.1  Background 

1.1.1  Linear  Ship-Motion  Theories 

During  the  last  twenty  years,  it  has  been  demonstrated  that  linear 
ship-motion  strip  theory  (for  example,  Salvesen,  Tuck,  and  Faltinsen,  1970) 
can  solve  with  good  accuracy  many  seakeeping  problems  related  to  the  average 
performance  of  ships  in  moderate  sea  conditions.  As  the  sea  condition  be¬ 
comes  more  and  more  severe,  however,  the  linearized  frequency-domain  strip- 
theory  approach  becomes  inadequate  for  analyzing  most  of  the  important  sea¬ 
keeping  problems.  More  importantly,  however,  the  linearized  approach  cannot 
be  used  to  forecast  some  of  the  most  important  responses  of  a  ship  even  in 
moderate  sea  conditions.  For  example,  relative  bow  motions  in  any  heading 
cannot  be  predicted  accurately  by  linear  strip-theory  ship-motions  computer 
programs . 


Strip  theory  not  only  linearizes  the  free-surface  boundary  condi¬ 
tions  and  the  ship-bull  boundary  condition,  but  it  also  replaces  the  three- 
dimensional  hydrodynamic  problem  by  a  summation  of  two-dimensional  sectional 
problems.  At  present,  it  is  not  clear  which  of  the  above  three  approxima¬ 
tions  of  linear  strip  theory  is  most  responsible  for  the  significant  errors 
in  the  predicted  ship-motion  results.  Ming  Chang  (1977)  has  developed  a 
fully  three-dimersional  ship-motion  theory  in  which,  similar  to  linear  strip 
theory,  both  the  free-surface  boundary  conditions  and  the  hull -boundary 
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conditions  are  linearized.  Some  of  the  unpublished  work  has  shown  that  there 
are  only  small  differences  between  the  heave  and  pitch  motions  predicted 
by  the  three-dimensional  theory  and  by  conventional  linear  strip  theory. 

These  preliminary  results  seem  to  indicate  that  the  inaccuracy  of  the  linear 
strip  theory  is  due  primarily  to  the  approximation  that  the  amplitude  of 
ship  motion  is  so  small  that  the  wetted  surface  of  the  hull  remains  fixed  over 
the  period  of  ship  motion.  Such  an  approximation  is  compatible  with  "wall- 
sided"  geometric  characteristics  of  ship  section  shapes,  even  for  moderate 
amplitude  ship  motions.  However,  most  ships  incorporate  a  considerable 
amount  of  flare  in  the  ship  sections  near  the  bow  and  stern.  Such  geometric 
characteristics  cause  substantial  nonlinear  hydrodynamic  effects,  particularly 
in  the  pitching  motions  of  the  ship.  These  nonlinear  effects  can  be  accounted 
for  by  modifying  the  linear  strip  theory  to  compute  the  hydrodynamic  character¬ 
istics  of  the  ship  as  a  function  of  time  and  allowing  for  the  instantaneous 
position  of  the  ship.  Furthermore,  the  linear  strip  theory  cannot  account 
for  the  changing  lateral  forces  and  moments  arising  from  the  changing  wetted 
surface  of  the  ship  as  it  heaves  and  pitches.  Hence,  the  important  coupling 
between  the  heave-pitch  and  the  sway-yaw  motions  is  missing,  but  this  coupling 
can  be  accounted  for  in  a  nonlinear  strip  theory  which  does  allow  for  changing 
wetted  surface  as  the  ship  oscillates. 

The  important  point  to  consider  is  that  the  motions  of  a  ship  are 
rarely  as  small  as  the  assumptions  of  linear  strip-theory  require.  For  example. 
Figure  1  below  shows  the  bow  motions  that  were  computed  using  linear  strip 
theory  for  a  destroyer  hull  in  regular  waves  with  the  very  modest  value  of 
steepness,  H/X  a  0.013  (where  H  is  the  wave  height  and  X  is  the  wave  length). 
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Figure  1:  Bow  Motion  of  Destroyer  Hull  in  Sinusoidal  Wave 

with  Wavelength,  X  =  1.20L,  Waveheight,  H/X  =  0.013, 
and  Froude  Number,  Fn  =  0.35  (Salvesen,  1978). 

Note  the  substantial  changes  of  bow  immersion  over  a  cycle  of  the  ship  motion. 
Clearly,  such  large  excursions  of  the  bow  relative  to  the  mean  level  of  the 
water  violates  the  basic  assumptions  of  linear  strip  theory.  It  is  intui¬ 
tively  obvious  that  such  large  changes  of  bow  immersion  will  entail  tremendous 
pitch-yaw  coupling  effects  that  are  completely  missing  in  the  standard  linear 
strip  theory  of  ship  motions.  It  is  precisely  this  aspect  of  the  linear 
strip  theory  that  we  feel  is  in  the  most  urgent  need  of  improvement. 

1.1.2  Nonlinear  Ship-Motion  Theories 

The  formulation  of  the  nonlinear  ship-mosion  problem  is  so  compli¬ 
cated  that  there  have  been  few  attempts  to  develop  computational  methods  for 
solving  such  fully  three-dimensional  problems.  Paulling's  (1974)  attempt  is 
one  of  the  few  exceptions.  He  assumed  that  the  frequency  of  encounter  is  low 
enough  that  all  he  needs  to  account  for  is  the  large  amplitude  hydrostatic 
effects.  He  uses  very  simple  approximations  to  the  hydrodynamic  effects. 

Paul  ling  has  had  some  success  with  this  method  in  determing  the  capsizing  of 
ships  due  to  low  cycle  roll  response  in  following  seas.  However,  this  sim¬ 
plified  method  cannot  be  used  to  predict  the  hydrodynamic  loading  on  ships 
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since  the  local  hydrodynamics  are  not  computed  with  the  accuracy  needed  to 
determine  local  pressure  distributions.  Salvesen  (1978)  has  formulated  a 
second-order  strip  theory  for  predicting  large  amplitude  ship  motions,  but  has 
not  presented  any  numerical  results  obtained  by  using  this  theory. 

There  is  a  need  for  a  new  computational  method  for  predicting  ship 
motions  which  takes  into  account  the  most  important  nonlinear  hydrodynamic 
effects  associated  with  large  amplitude  motions.  However,  little  progress 
has  been  made  towards  this  goal  since  it  is  so  difficult  to  reduce  the  compli¬ 
cated  nonlinear  three-dimensional  mathematical  formulation  to  a  simpler  form 
from  which  a  computationally  feasible  computer  program  can  be  developed. 

This  requires  the  determination  of  those  nonlinear  effects  which  must  be  in¬ 
cluded  and  those  which  can  be  disregarded  without  a  substantial  loss  of  the 
accuracy  required  to  predict  the  observed  ship  motions. 

1 .2  Summary  of  New  Nonlinear  Method 

1.2.1  Objective 

The  objective  of  this  work  has  been  to  develop  a  new  method  and  appro¬ 
priate  computer  program  for  predicting  the  nonlinear  large-amplitude  motions 
for  a  "non-wall -sided”  ship  advancing  in  a  seaway.  The  local  hull  form  above 
and  below  the  still  water  has  been  accounted  for  so  that  phenomena  such  as 

bow-flare  impact,  slamming  and  the  shipping  of  water  on  deck  as  well  as  the 

ship  motions  themselves  can  be  predicted  with  fair  accuracy  in  relatively 
severe  sea  conditions. 
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1.2.2 


Theoretical  Method 


A  new  nonlinear  time-domain  "strip-theory"  calculation  method  for 
ship  motions  is  outlined  in  Section  2.  In  this  method  it  is  assumed  that  the 
frequency  of  encounter  is  sufficiently  small  so  that  the  inertia,  damping, 
and  diffraction  forces  are  much  smaller  than  the  hydrostatic  restoring  and 
the  Froude-Kriloff  exciting  forces.  The  added-mass,  damping,  and  diffraction 
terms  in  the  equation  of  motion  are  then  predicted  by  the  conventional  linear 
ship-motion  theory,  whereas  the  hydrostatic  restoring  and  Froude-Kriloff 
forces  are  computed  by  integrating  the  pressure  over  the  instantaneous  wetted 
hull  surface  at  each  time  step. 

The  low-frequency  assumption  restricts  the  applicability  of  the 
present  method;  however,  it  should  be  suitable  for  solving  beam,  quartering, 
and  following  sea  problems.  It  should  also  be  applicable  to  bow  and  head 
sea  cases  with  low  forward  speed. 

The  most  important  reason  for  applying  the  low-frequency  assumption 
and  for  using  the  linear-theory  added-mass,  damping,  and  diffraction  terms 
is  that  this  is  a  natural  first  step  towards  a  more  complicated  nonlinear 
theory  where  the  total  two-dimensional  sectional  hydrodynamic  problem  is 
solved  at  each  time  step.  Furthermore,  as  shown  in  the  results,  the  low- 
frequency  approach  gives  us  some  good  realistic  estimates  of  the  magnitude  of 
some  of  the  nonlinear  effects  not  included  in  the  conventional  linear  strip- 
theory  approach. 
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Section  2 


THEORETICAL  APPROACH 

A  ship  advancing  at  constant  mean  forward  speed  in  oblique  waves 
is  being  considered.  It  is  assumed  that  the  frequency  of  encounter  is  suf¬ 
ficiently  small  so  that  large  amplitude  body  motions  will  result  only  in 
small  free  surface  disturbances.  Let  ( x,y,z )  be  a  right-handed  orthogonal 
coordinate  system  fixed  with  respect  to  the  mean  position  of  the  body,  with 
z  vertically  upward  through  the  center  of  gravity  of  the  body,  x  in  the  direc¬ 
tion  of  forward  motion,  and  the  origin  in  the  plane  of  the  undisturbed  free 
surface. 


2.1  Incident  Wave 


The  incident  wave  is  represented  by  a  second-order  Stokes  wave  with 
the  velocity  potential  given  by 


a  -  *0)  +  *(2)  +  -  iga.  a-ikX  kz  iwt 

$1  !  ?]  +  9 j  +  . . .  =  e  e  e 


(2.1) 


and  the  wave  elevation  by 


5(X,t)  =  a  cos(kX  +  wt)  +  hka^  cos2(kX  +  wt)  +  ...  (2.2) 


where 


X  =  x  cosg  +  y  sing. 


Here  a  is  the  first-order  wave  amplitude,  a  is  the  wave  frequency,  k  is  the 
wave  number,  g  is  the  heading  angle,  and  co  is  the  frequency  of  encounter. 
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Note  that  the  expression  (2.1)  is  correct  to  the 

wave  amplitude,  and  we  are  using  a  notation  such  that 
(2)2 

'  is  of  0(a  )  where  a  is  a  small  parameter  proportional 


second-order  in 
is  of  0(a)  and 
to  the  wave  slope. 


2.2  Velocity  Potential 

It  is  assumed  that  the  total  velocity  potential,  $(x,y,z,t)  has 
the  following  expansion 

$  =  $0)  +  +  ...  (2.3) 

where  the  body  disturbance  potential,  is  of  0(a5)  where  <5  is  a  small 

parameter  related  to  the  low  encounter  frequency  and/or  the  slenderness  of 
the  body. 


2.3  Equations  of  Motions 

The  equations  of  motion  may  be  expressed  as 


6 


2 


-M . .  ri. 

Ji  1 


r.?0J  *  ,  j 

J  J 


(2.4) 


where  i  or  j  equal  to  1 ,  2,  3,  4,  5,  6  refer  to  surge,  sway,  heave,  roll, 
pitch,  and  yaw,  respectively,  and  where  Mjj  is  the  generalized  mass  matrix, 
n.j  are  the  components  of  the  acceleration  of  the  center  of  gravity.  Further- 
more,  Fj  are  the  forces  and  moment  components  resulting  from  the  velocity 

.  .  add 

potential  (2.3),  and  F.  are  additional  force  and  moment  components  due  to 

J 

viscous  drag,  lift,  and  vortex  shedding  effects.  Note  that  the  moments  are 
about  the  axis  of  a  coordinate  system  parallel  to  the  (x,y,z)  system,  with 
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the  origin  at  the  center  of  gravity  of  the  body. 


2.4  Hydrodynamic  Forces 


The  hydrodynamic  force  and  moment  components  on  the  body  are 

f/°T  =  "st  pnj  ds  (2-5) 

where  n.  are  the  components  of  the  generalized  normal  and  St  is  the  instan- 
j  1 

taneous  wetted  hull  surface.  The  pressure,  p,  is  by  Bernoulli's  equation 


p  =  -pH "  I  " pgz 


(2.6) 


If  we  only  include  terms  of  0(a),  0(a),  and  O(afi),  the  second-order  hydro- 
dynamic  force  due  to  the  velocity  potential  (2.6)  alone  is 


F/0T  =  (p^L*j  +  pgz)  n.  ds  (2.7) 

p(»s|V*J1)|2  +  ^‘41,1))  n.  ds 

Here  the  first  integral  is  over  the  instantaneous  wetted  surface,  S whereas 
the  last  integral  is  over  the  undisturbed  position,  SQ. 

POT 

It  is  convenient  to  decompose  the  force,  F.  into  the  following 

J 

components 


F.TOT.FiFK  +  F1R  +  F.m  +  F.D*F.DC 
0  J  J  J  J  J 


(2.8) 


Fk 

where  by  equation  (2.6)  the  Froude-Kriloff  exciting  force,  Fj  and  the  hydro- 

R 

static  restoring  force,  Fj  are 
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(2.9) 


FK  _  R 


Fj  +  Fj  =  -ffs  4 1  +  pgz^  nj  ds  * 

and  the  hydrodynamic  forces,  F.^,  and  the  diffraction  exciting  force, 

J  J 


are 


p  HD  ,  c  D  _  ,,  3  .O-1)  „  . 

Fj  +  Fj  -  "St  *8  "j  ds  • 

and,  finally,  the  second-order  DC  force  is 


(2.10) 


Fj0C  ■  -"s0  f  lV1o')|2  "J  ds  • 

Since  the  body  disturbance  potential,  is  of  first  order  in 

a,  this  potential  can  be  assumed  to  be  harmonic.  Hence,  we  may  write 


(2.11) 


♦J1 •' 1  ■  <«j+J  *  ♦D>e'“t 


(2.12) 


where  4>q  is  the  diffraction  potential,  and  <j>j  is  the  contribution  to  the  poten¬ 
tial  from  the  jth  mode  of  motion  By  substituting  (2.11)  into  (2.9)  it 
follows  that 


F. 

J 


HD 


n.  +  F . 


(2.13) 


where  A.,  and  B , ,  are  the  added  mass  and  damping  coefficients.  Within  the 

J  I  J  * 

present  formulation  of  the  problem,  the  terms,  Aj^,  Bjj,  and  FjD,  can  be 
computed  by  linear  ship  motion  theory.  Note  that  if  strip  theory  gives  in¬ 
adequate  estimates  for  these  terms  at  low  frequencies,  a  three-dimensional 
method  may  be  used  for  the  prediction  of  these  linear  hydrodynamic  terms. 


The  total  hydrodynamic  force  on  the  body  consists  of  the  parts  of 

the  potential  force,  F.^,  and  the  additional  force,  F-^,  which  has  the 

J  J 
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following  components 


ADD  .  f  VD  +  L  +  K 
J  J  D  0 

Here  F.v^,  F and  F.  ^  are  the  forces  due  to  viscous 
3  J  j 

shedding,  respectively. 


drag. 


(2.14) 
lift,  and  vortex 


2.5  Problem  Formulation  Summary 

The  equations  of  motion  which  are  to  be  solved  by  a  finite-dif¬ 
ference  time-simulation  method  are 


6 

I 

i=l 


V  =  [F«  +  F.R]  +  (FU  -  -  B..n.] 

Zi  Ji  i  j  j  j  ji  i  ji  i 


(2.15) 
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3  3  3  3 


FK  R 

where  F.  +  F.  are  the  nonlinear  Froude-Kriloff  and  restoring  forces 
J  3 

given  by  equation  (2.8),  [FjD  -  A^fy  -  are  the  linear  diffraction, 

added  mass,  and  damping  terms,  Fj  D  is  the  second-order  DC  force  given  by 
equation  (2.11 ),  and  FivD>  p .  ,  and  F.  are  the  viscous  drag,  lift,  and 

J  J  J 

vortex  shedding  terms. 
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Section  3 


COMPUTATIONAL  METHOD 

3.1  Completed  Computer  Codes 

Two  computer  codes  have  been  developed  which  determine  the  non¬ 
linear  large-amplitude  ship  motions  in  the  time  domain.  One  computer  code 
predicts  the  heave,  sway,  and  roll  motions  for  a  ship  in  beam  seas,  whereas 
the  other  code  predicts  the  heave  and  pitch  motions  in  head  and  following 
waves.  In  both  of  these  codes,  the  incident  wave  is  a  second-order  Stokes 
wave  (see  equations  (2.1)  and  (2.2)).  The  diffraction  force,  F^,  and  the 
added-mass  and  damping  coefficients.  A- .  and  B. . ,  respectively,  are  predicted 

'  J  '  J 

for  the  particular  frequency  of  encounter  (note  that  the  nonlinear  Stokes 

wave  is  periodic)  by  linear  ship-motion  strip  theory  of  Salvesen,  Tuck,  and 

Faltinsen  (1970).  The  viscous  damping  terms  are  only  included  by  an  approxi- 

FK 

mate  value  for  roll  damping.  The  Froude-Kriloff  force,  F.  ,  and  the  dif- 

J 

fraction  force,  FjD,  are  computed  at  each  time  step  by  integrating  the  pres¬ 
sure  over  the  instantaneous  wetted  surface  (see  equation  (2.9)). 

3.2  Time  Stepping  Procedure 

The  motions  of  the  ship  are  computed  as  a  function  of  time  by  a 
numerical  procedure  in  which  the  time,  t,  is  advanced  in  small  time  steps. 

At.  We  shall  suppose  that  the  displacements,  velocity,  and  acceleration  com¬ 
ponents,  nj(t),  >ij(t),  and  fi j ( t ) ,  respectively,  are  known  at  time,  t.  When 
the  time  is  advanced  from  t  to  t  +  At,  the  new  displacement  and  velocity 
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components,  can  be  determined  from  the  known  motions  as  follows: 

n j  (t+At)  =  nj(t)  +  fij(t)  •  At  +  ^(t)  •  (At)2  (3.1) 

and 

iij  (t+At)  =  rij(t)  +  Tij(t)  •  At  .  (3.2) 

Knowledge  of  the  displacements,  rij(t+At)  and  the  velocities, 

J 

n.(t+At)  allows  one  to  compute  the  force  component  on  the  right-hand  side  of 

J 

the  equation  of  motion  (2.15).  One  can  then  use  the  equation  of  motion  (2.12) 
to  compute  the  updated  value  of  acceleration,  n^t+At). 

%J 

This  time  stepping  procedure  seems  to  work  well  for  the  beam  sea 
cases  investigated.  However,  we  encountered  instabilities  with  this  approach 
for  following  sea  cases,  and  then  introduced  a  second-order  method  where  the 
displacement  at  time,  t+At,  is  given  by 

Hj(t+At)  =  nj(t)  +  f|j(t)  At  +  1  [2rij (t)  +  ri j ( t+At )]  (At)2  (3.3) 

and  the  velocity  by 

n  ( t+At)  =  f|.(t)  +  H  [ri(t)  +  fi( t+At)]  At  .  (3.4) 

J  J 

Note  that  the  displacements  and  the  velocities  at  time,  t  +  At,  are  here  ex¬ 
pressed  in  terms  of  the  unknown  acceleration,  n(t+At),  so  that  in  this  case, 
an  iteration  scheme  has  to  be  used  for  each  time  step.  This  approach  resulted 
in  stable  solutions  for  most  of  the  cases  investigated  with  At  equal  to  one 
twentieth  of  the  period.  However,  we  found  that  for  some  particular  following 
sea  cases  with  extremely  low  encounter  frequency,  it  was  difficult  to  reach 


stable  conditions.  Other  time  stepping  procedures  will  be  considered  in 
future  work. 
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Section  4 


RESULTS 


Two  ship  hull  forms  were  used  to  evaluate  and  to  explore  the  poten¬ 
tials  of  the  two  new  time-domain  nonlinear  large-amplitude  computer  codes  for 
predicting: 

(i)  heave,  sway,  and  roll  motions  in  beam  seas;  and 

(ii)  heave  and  pitch  motions  in  head  and  following  seas. 

The  two  hull  forms  used  were  the  Series  60  form  with  block  coefficient, 

Cg  =  0.70,  and  a  typical  trawler  form.  The  body  plans  for  the  Series  60 
and  the  trawler  form  are  shown  in  Figures  2  and  3,  respectively.  Eleven 
stations  were  used  for  each  ship  form.  The  selected  segments  used  for  each 
station  are  shown  in  the  figures. 

These  two  hull  forms  were  selected  since  the  Series  60  form  is 
almost  "wall -sided",  whereas  the  trawler  form  has  extreme  "V"  sections.  It 
was  felt  that  by  using  these  hull  forms  the  nonlinear  effects  strictly  due  to 
the  large-amplitude  motions  could  be  separated  from  the  effects  strictly  due 
to  the  "non-wall -sidedness". 

4.1  Beam  Sea  Results 

The  accuracy  of  the  beam  sea  computer  code  was  first  evaluated  by 
using  small  amplitude  waves  (wave  amplitude/draft  ratio,  a/D  *  0.01).  The 
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Fioure  3:  Body  Plan  of  Trawler  Form  with  the 
Selected  Serpents  for  each  of  the 
Eleven  Stations. 


pitch,  heave,  and  roll  motions  predicted  by  the  time-domain  method  were  found 
for  most  of  the  wavelength  cases  investigated  to  agree  within  1%  of  the  results 
obtained  by  conventional  linear  frequency  domain  strip  theory.  For  larger 
wave  amplitudes,  we  found  discrepancies  between  the  nonlinear  time-domain  code 
and  the  linear  frequency  domain  theory.  Som  sample  results  for  the  largest 
wave  amplitude  investigated  (a/D  =  0.50)  are  shown  in  Figure  4.  For  this 
case,  the  maximum  roll  displacement  is  29.7  degrees,  and,  as  seen  in  the  fig¬ 
ure,  the  deck  is  partly  submerged.  The  computed  time-domain  motions  were 
periodic  and  stable  for  this  large  amplitude  case,  and  the  motions  were  about 
20%  different  from  the  linear  theory  results.  Computations  were  also  made 
with  larger  wave  amplitude;  however,  we  then  found  that  the  results  were 
somewhat  unstable.  These  effects  need  to  be  further  investigated. 

Figure  5  shows  the  sway  displacement,  n2  computed  by  the  time-domain 
code  as  a  function  of  time  for  the  Series  60  hull  in  beam  waves  with  a/D  =  0.25. 
The  results  clearly  show  the  nonlinear  drift  which  cannot  be  predicted  by 
linear  theory.  The  sway-displacement  response  shown  in  Figure  5  is  approxi¬ 
mately 


n2  =  A,]t  +  A2  cos 

where  A^t  is  the  nonlinear  drift,  and  Ag  cos  wt  is  the  harmonic  motion.  The 
linear  theory  predicts  a  sway  amplitude  of  0.2035  for  this  case,  whereas  for 
the  nonlinear  time-domain  results  shown  in  Figure  5,  we  have  that  0.2040  < 

A2  s  0.2055.  These  seem  to  indicate  stable  results  which  are  in  good  agree¬ 
ment  with  linear  theory. 
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Figure  4:  Comnuted  Large-Amplitude 
Roll  Motion  in  Beam  Seas 
for  Series  60  Hull. 


Sway  Displacement, 
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Head  and  Following  Sea  Results 


The  heave  and  pitch  responses  in  head  and  following  seas  were  in¬ 
vestigated  for  both  the  Series  60  hull  and  the  trawler  form  for  several  com¬ 
binations  of  wave  conditions  and  forward  speeds.  For  small  wave  amplitudes 
(a/0  =  0.01),  we  found  good  agreement  with  linear  theory  and  stable  results 
other  than  for  following  sea  cases  with  the  frequency  of  encounter  close  to 
zero.  For  the  very  low  frequency  cases,  the  results  were  clearly  unstable. 

i.ci -ge  amplitude  responses  were  also  investigated  and  Figure  6  shows 
the  difference  in  the  heave  displacements  between  the  conventional  linear 
strip  theory  and  our  new  nonlinear  time-domain  computer  code  for  the  trawler 
at  zero  forward  speed  in  head  waves  (ship  length/wavelength  ratio,  L/X  =  0.75). 
It  is  seen  that  for  a  moderately  steep  wave  with  waveheight-to-wavelength  ratio, 
H/X  =  0.03,  the  difference  is  15%  for  the  downward  heave  displacement  and  5% 
for  the  upward  displacement.  For  a  steep  wave  with  H/X  =  0.09,  the  differ¬ 
ences  are  36%  and  19%,  respectively. 

In  Figure  7,  the  corresponding  results  are  presented  for  the  Series 
60  hull.  The  abrupt  changes  in  the  two  curves  are  due  to  bow  submergence  at 
H/X*  0.035.  The  effect  of  the  water  on  deck  is  included  in  the  computer  code; 
however,  it  is  questionable  if  the  method  in  its  present  form  gives  a  meaning¬ 
ful  representation  of  the  actual  hydrodynamic  effects  associated  with  the 
large  amount  of  water  passing  over  the  deck.  This  problem  area  needs  to  be 
further  investigated. 
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Figure  7:  Differences  in  Heave  Displacements  Between 
Linear  and  Nonlinear  Theory  for  Series  60, 
CQ  =  0.70,  at  Zero  Forward  Sneed  (L/X  =  0. 


Waveheight-to-Waveleng^h  Ratio,  H/X 


Section  5 


CONCLUDING  REMARKS 


The  results  obtained  by  the  two  computer  codes  for  predicting  non¬ 
linear  large-amplitude  motions  for  ships  in  beam  seas  and  in  head  and  fol¬ 
lowing  seas  have  been  found  in  the  linear  range  to  agree  well  with  linear 
frequency-domain  strip  theory.  For  steep  waves,  it  has  been  found  that  there 
are  large  differences  between  the  nonlinear  method  and  the  linear  theory. 
Further  investigations  are  needed  in  order  to  determine  the  accuracy  of  the 
nonlinear  effects  determined  by  this  method,  but  the  results  clearly  demon¬ 
strate  that  nonlinear  effects  are  of  significant  magnitude  for  large  amplitude 
motions. 


It  is 

can  be  achieved 
and  diffraction 


strongly  believed  that  improved  accuracy  in  the  nonlinear  range 
by  including  the  nonlinear  effects  for  the  added-mass,  damping, 
terms  in  addition  to  the  restoring  and  Froude-Kriloff  forces. 
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